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Abstract 

Multisite interaction spin-S models in an external magnetic field are studied recursively 
on the Bethe-like lattices. The transfer-matrix method is extended to calculate exactly 
the two-spin correlation functions. The exact expressions for the correlation length and 
magnetic susceptibility are derived for spin-1/2 models. The singularity of the correlation 
length with critical index v = 1 and the proportionality of magnetic susceptibility to 
correlation length in the second order phase transition region of spin-1/2 ferromagnetic 
models on the Bethe-like lattices are established analytically. 

1 Introduction 

The investigation of physical systems such as binary alloys U, classical fluids |§, liquid 
bilayers [|J, solid 3 He ||, rare gases |5j and anisotropic magnets such as CeBi, EuSe, etc. lead 
to multisite interaction (MSI) spin and Heisenberg models. This models have a complex phase 
diagrams and one of the challenging problems in modern statistical mechanics is the effect 
of multisite interactions on the properties of order-disorder and order-order phase transitions. 
Mean-field theory which is most used and the simplest method for obtaining approximate re- 
sults, has been shown to give poor results for such systems H |6], 0. Therefore, to obtain better 
approximations high temperature series expansions ||, Monte Carlo simulations 0, Bethe lat- 
tice and Cayley tree approximations |8[]-|[10|| are used. Although, there are lots of approximate 
methods one need to have simple but accurate technique to approximate multisite interac- 



tion (MSI) spin systems. After the pioneer works of Bethe and Peierls |TT| the Bethe lattice 
approximation has attracted particular interest because the calculations on the Bethe lattice 
are rather simple and results reflects essential features of the systems, even when conventional 
mean- field theories fail 0]. The main difference between the Cayley tree and Bethe lattice 
as discussed in Ref. [13] is that in Cayley tree the surface plays a very important role. The 
number of sites on it's surface is a finite fraction of the total sites even in the thermodynamic 
limit. As a consequence, the models on the Cayley tree exhibit quite unusual critical behavior 



without long-range order P]-[l7[]. To overcome the problems with boundary one considers only 
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properties of sites deep within the interior of the Cayley tree. The union of such equivalent 
sites is regarded as forming Bethe lattice and it is assumed to have a translational symmetry 
like any regular lattice. 

It was shown by Monroe H that one can obtain much better approximations than on the 
Bethe lattice by using Husimi or more general Bethe-like lattices to investigate MSI systems, 
especially in case of frustration. By going to Bethe-like lattice one constructs the system from 
different type of blocks and, hence, can differentiate between two lattices having the same 
coordination numbers, such as triangle and simple cubic lattices, whereas in the case of Bethe 
lattice it is impossible. Furthermore, the various options regarding the connection of the basic 
building blocks, the symmetry of Hamiltonian and, as a consequence, the possibility of division 
of such lattices on shells are of importance (see Sect. 2). 

In the investigation of critical phenomena the correlation functions play an important role. It 
is established that near second order phase transition critical temperature T c the fluctuations in 
system rise and grappled all the system at T c . It is believed that singularities of thermodynamic 
quantities such as the specific heat and the susceptibility at T c are related to the increase of 



fluctuations in the system or, in other words, to divergence of correlation length at T c [jTS 



Since two-spin correlation functions can be measured experimentally as linear response to an 



adiabatic or isotermic applied field, or scattering of neutrons or electromagnetic waves ||19||-[ 
it is desirable to have an analytic form of the correlation function to locate phase transitions 
and a deeper understanding of phase transition phenomena. 

Recently spin-spin correlation functions of spin-5* models in an external magnetic field was 



calculated on the Bethe lattice [[21J , where the crucial role of Bethe lattice dimensionality was 
demonstrated in determining the critical behavior of correlation length. It was shown that 
correlation length diverges at the critical point of m and \ for the Bethe lattice Ising model 
with critical exponents v — 1, which is different from the mean- field critical exponent v — |. 

In this paper we apply the transfer matrix method to analytical calculation of the two- 
spin correlation functions of the multisite interaction spin-S 1 models in the presence of an 
external magnetic field on the general Bethe-like lattices with Hamiltonian symmetric under 
the permutation of spins at any p-polygonP] of the lattice. The derivation of the analytic 
expression for correlation functions for arbitrary spin-S* is presented in Sect. 2. For spin-| case, 
taking into account the dimensionality of the Bethe-like lattice, the exact expressions of the 
correlation length £ and magnetic susceptibility \ are derived in Sect. 3. It is demonstrated 
exactly that at T = T c the correlation length £ diverges with critical exponent v — 1 and x ~ £ 
for arbitrary p-polygon Bethe-like lattice with q p-polygons going out from each site. 



2 The Model 

The Bethe-like lattice may be constructed as follows : Connecting (q — 1) p-polygons at a 
single base site one gets a first generation branch (Fig.l.(a)). Connecting (q— 1) first generation 
branches to each site of a new p-polygon (basic block), except the base site of this system 
one gets the second generation branch (Fig. l.(b)). Continuing this process one gets higher 
generation branches. To complete the system, one connects q nth generation branches at their 

V-polygon is a closed polygon with p edges (sites), for example 2-polygon is a line, 3-polygon is a triangle 
etc. 
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base site (Oth or central site) (Fig. l.(c)). The set of sites deep within this system is called the 
Bethe-like lattice (see the Introduction). As we see only topological and connectivity properties 
are used to define the Bethe-like lattice. This system is not a lattice in a conventional sense. 




Figure 1: The recursive structure of 4-polygon Bethe-like lattice with q = 3. (a) The first 
generation branch; (b) the second generation branch; (c) the central part of 4-polygon Bethe- 
like lattice with q = 3, the numbers 0,1,2 indicate shells. 

The Bethe and Bethe-like lattices are usually viewed as ramified trees embedded in infinite 
dimensional Euclidean space, with a constant vertex connectivity. The fact that Bethe and 
Bethe-like lattices usually are viewed as "infinite-dimensional" follows from the definition of 
lattice dimensionality d in Euclidean space given in Ref. [13] 

d = lim - — -, (1) 
where C n is the total number of sites within n steps of a given site. For a Bethe-like lattice 



C n = l+Y,C m , C n = q(p-l)[(q-l)(p-l)} n -\ C„ = l, (2) 

m=l 



where C n is the number of sites on the nth shell. The substitution of this expression into ([[]) 
gives d = oo, so in this sense the Bethe-like lattice is "infinite-dimensional". 

Recently it was shown that Bethe-like lattices can be embedded in two-dimensional space 
of constant negative curvature (the H2 hyperbolic plane) with fixed bond angles and lengths. 
This hierarchical structures like Bethe-like lattices can be described as a conformal tiling of the 



hyperbolic plane and their metric properties are studied |22 . 

After this brief survey of the geometrical properties of Bethe-like lattices let us now define 
the Hamiltonian of the multisite interaction (MSI) system. The general MSI spin Hamiltonian 
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have the form 

n = - J f p E IlS;-§£^ (3) 

(polygons) i 

where Si takes values S, S — 1, . . . , —S + 1, — S; the first sum goes over all the p-polygons of the 
lattice and II means the product of all spins placed on a p-polygon; the second term presents 
the interaction with an external magnetic field and the sum goes over all the sites of the lattice. 

The advantage of the Bethe-like lattice is that for the models formulated on it, the exact 
recurrence relation can be derived and the dynamical systems theory may be used to explore 
the thermodynamic properties of the models || [10], |23[. Let us denote the partition function 
of the N — th generation branch with the basic spin in state S as B(N, S). By cutting the 
lattice at the basic p-polygon one will obtain (q — l)(p— 1) interacting (N — 1) — th generation 
branches 

B(N, S ) = £ w(S , S 1 ,..., S^B^N - 1,S 1 ) ■ ■ ■ B^(N - 1, Vi), (4) 

S\,...,S p -i 

where So, Si, ... , S p -i are the spins at the basic p-polygon, w(So, Si, ... , S p -i) is the statistical 
weight of the basic p-polygon and 7 = q — 1. 

w(S , Si, ... , S p _i) = exp |^S Si ■ ■ • S p _i + ^(S H h S p _i)| , (5) 

where (3 = 1/kT. By cutting apart the lattice at the central site one can obtain for the partition 
function Z = e _/3W the following expression 

Z N = Y,B q (N,S ). (6) 

So 

From Eq. @ follows that w(Sq, Si, ... , S p _i) is a symmetric function under the permutation of 
its arguments So, Si, ... , S p _i. Now it is easy to check that all the spins at the n — th generation 
branch surface are equivalent and we refer to them as n — th shell (see Fig. l.(c)). 
Let us now define the auxiliary quantities x(s) as follows 

xn{s) = ^7 S J. , x N (S) = l and s = S - 1, . . . , -S. (7) 
ti(l\,b) 

The thermodynamic quantities such as the magnetization, the specific heat etc. and even 
thermodynamic potentials such as the free energy p4| may be expressed in terms of x(s). Of 
course, x(s) has no direct physical meaning, but {x(s)} determine the state of the system in 
the thermodynamic limit. From Eqs. (Q) and (0) we can get 2S recurrence equations for 

Xn = Eft,...^.! Si,..., S p _ 1 )x^(S 1 ) ■ ■ ■ x^(S p _i) 
N+1 Esu^Sp-! W (S, Si, ... , S p _i)x^(Si) ■ ■ ■ a;^(Sp_i) ' 

Using Eqs. @ and (0) one can express the magnetization of the central site as follows 

»~\\S) I' S» Zs x%(S ) ' W 
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where (...) denotes the thermal average and \i takes values from 1 to 25'. 

Further we will restrict our treatment to the case when the series of solutions of recurrence 
equations (§) converges to a stable fixed point as N — > oo. Note that for stable fixed point of 
Eq. (H) the difference between xjy(s) and xn-i{s) disappears for N — > oo 

lim XN-i(s) = x(s) 

for all finite values i. The corresponding fixed points equations of the mapping (§) have the 
form 

[w(s, Si,..., Sp-i) - x(s)w(S, Si,..., Sp-i)] x 7 (S'i) • ■ ■ x 7 (S p -i) = 0, (10) 

Si 

where % = 1, . . . ,p — 1. 

The correlation functions are defined as follows 

g(n) = ^((S S n )-(S }(S n }), (11) 

where spin So belongs to the central polygon and S n located somewhere in the nth shell. Using 
denotation (^) and recurrence relations (|) we can express (SoS n ) as follows 

(S S ) = lim ^ S ° x ~h( S o) w ( S Qiyii ■ ■ ■ z i) x N-i(yi) ■ ■■ x N-i( t i) x 
°™ N ^°° E^{Si)w(Si,yi,...,ti,zi)xJ f _ 1 (yi)---x 7 N _ 1 (ti)x 

x]i\(zi)w(zi } y 2 , . . . ,t 2 , z 2 )x' y N _ 2 (y 2 ) ■ ■ ■ xlf_ 2 (t 2 ) x ■■■ 

x J N ~\(z 1 )w(zi,y 2 , ...,t 2 , z 2 )x^_ 2 (y 2 ) ■ ■ ■ x J N _ 2 (t 2 ) X ••• 

xJi_ n+1 (z n -i)w(z n -i, y n , . . . , t n , S n )xjf_ n (y n ) • • • xJf_ n (t n )xJj^ n (S n )S n 

xjj- n +l( z n-l) w ( z n~l, Vni ■ ■ ■ i^ni S n )xjf_ n (y n ) ' ' • xJf_ n (t n )xJf_ n (S n ) 

where sums run over all spin states displayed. 

To use the transfer matrix method let us introduce the following symmetric matrixes with 
elements 

7-1 

Mi{zi,z i+ i) = [x N -i{zi)x N ^i{zi + i)\ 2 x 

53 w ( z i,yi+u---,U+i,z i+1 )xJj_ i _ 1 (y i+1 )---x' y N _ i _ 1 (t i+1 ) 1 (13) 

A(S n ,S ) = [x^niSnMSo)]^ , (14) 

7+1 

A(S n ,S ) = S n S [x N „ n (S n )x N (S )} 2 . (15) 

For stable fixed points {x(s)} of Eq. (§D all matrixes Mi coincide, Mi = M for all i. Thus, we 
can rewrite {SoS n ) as follows 

{SoSn) - Tr(AM n ) ■ (16) 

Let us now briefly discuss some properties of Mi matrixes. Above we stated that for stable fixed 
points {x(s)} all the matrixes M, coincide, but this is true only when w(Sq, Si, ... , S p -i) is a 
symmetric function under the permutations of its arguments, which holds for MSI Hamiltonian 
(|3|). In the case when one wants to introduce also pair interactions one find that the definition 



(12) 
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of Mj becomes not unique. To demonstrate this let us consider spin-| Ising model on the 
4-polygon (square) Bethe-like lattice with J2 pair interactions only. In Fig. 2 one can see 
two configurations of spins on a 4-polygon which differ only in permutation of spins, but their 
statistical weights are different. For example, for configuration (1) we have 

M (+,+) = w(+,y i+1 ,t i+1 ,+)x' y (y i+1 )x 7 (t i+1 ) = 

and for configuration (2) respectively 

M(+, +) = e 4/3Ja + e" 4/3j2 x 27 + 2x\ 



+ 



+ 




+ 




w(+, -,-,+) = ! w( 



) ) ) 



+ = e 



-4(3J 2 



Figure 2: Two configuration of spins on a 4-polygons with asymmetric statistical weight func- 
tion. 

To recover the symmetry of w we can add all possible diagonal interactions between the 
spins on a polygon (see Fig. 3). Thus, if one wants to consider pair or /and /-interactions 
(/ = 2, . . . ,p — 1 is the number of interacting spins on a p-polygon) one has to consider all 
possible interactions between spins on a polygon to have a symmetric function of the statistical 
weight w under the permutation of spins on a p-polygon. 

In order to use the transfer-matrix method we consider only the case when w is a symmetric 
function as described above. Now let us return to matrices M and A. Using the symmetry 
property of w, Eq. (f|) and definitions flllf ), ([14]) one can easily prove that matrices A and M 
commute. Thus they can be diagonalized simultaneously. The spectrum of matrix A is highly 
degenerated, it has only one nonzero eigenvalue 

s=S-l,...,-S 

with eigenvector 

After the diagonalization of matrix M we come to the following expression for the correlation 
functions 

g(n) = T [ A k [^) , (17) 
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Figure 3: The schematic of the 4-polygon Bethe-like lattice with all possible pair interactions 
on a polygon. 



where A& are the eigenvalues of matrix M. The coefficients Ak (k = 1, 2, . . . , 2S) can be obtained 
by solving the following system of linear algebraic equations 

{ g(0) = A 1 +A 2 + --- + A 2S , 
0(1) = Aih + A 2 l 2 + ■■■ + A 2S l 2S , 

{ 9(2S - 1) = All 3 - 1 + A 2 lf~^ + ■■■ + A 2S l? s -\ 
where = Xk+i/Xi- The solution of this system of equations can be presented in the form 

A ^ = — ^2sTj 7n > » k- (18) 

lli=il'* - k) 

where Fj(k) are the elementary functions of 25 — 1 variables li,l 2 , . . . , h-i, h+i, ■ ■ ■ , hs- 

F 2S -x{k) = 1, 

F 2 s-2(k) =h-\ lk-i + h+i H 1- ^25, 



F (k) — l\ • • ■ h-ih+i 4 • - hs 

(for more details see the Appendix of Ref. [21]). From the Eq. fll8D we can find for MSI spin-| 
Ising model (Sq) = ~ and 

A, = g(0) = 4 (S 2 ) - ml = 1 - m?, (19) 
for MSI spin-1 Ising model respectively 

Z 2 £(0) - 0(1) A 3 (m 2 - m?) - A 1 ((S S 1 ) - m?) 



A3 — X 2 



and A 2 = A\(lx 12), where m 2 = (So), etc. 

Thus, Eqs. (PD,(|TDD and flTT| ) give us a full set of equations for investigation of spin-S* MSI 
and pair interaction models on the Bethe-like lattices. 



3 The General MSI Spin-± Ising Model 

Here we apply the general technique developed in previous section to investigate the MSI 
spin- 1 Ising model in an external magnetic field h with Hamiltonian 



H = -2 P J P E i[s p -u 2 E S i S j -2h y £S i 

<polygons> <ij> i 



(20) 



where Si takes values ±|, the first term describes the multisite interactions on p-polygons, the 
second sum goes over all possible pairs of spins on each p-polygon (see Sect. 2) and the third 
one goes over all spins on a lattice. The statistical weight of a p-polygon has the form 



w(S , . . . , Sp-i) = ex P 



<ij>eP Q ieP 



(21) 



We see that w is a symmetric function under the permutation of its arguments. Let us make 
the following notations 



a k 




(22) 



2 q 



where k is the number of spins on a p-polygon in the state +|, k < p. Using this notations and 
xn(— ) = xn we can present the corresponding recurrence equations (Bj) as follows 



(23) 



where C p = p\/ {k\{p — k)\) is the number of ways of choosing k elements from p ones without 
permutations. 

The stable fixed points x of the mapping can be obtained from the fixed point equation 



and the stability condition 



p-i 

E 

fe=0 



(24) 



d_ 

dx 



f{x)\x=x c 



< 1. 



Using the identity C k p _ x = C k p _ 2 + C k Z 2 one can rewrite the fixed point equation (pi[) as follows 

p-2 



E (^ 7+1 ^+i - ^ + ^«fc+ 2 - ajk+i) c*_ 2 a^- 2 - fe ^ = 0. 

fe=0 



(25) 



In the stable fixed point the magnetization of the central site has the form 

l-x q 

mi = . 26 

1 + xi y ' 

It is well known that a ferromagnetic system at high temperatures (T > T c ) is in the paramag- 
netic phase with m x = and, hence, x=l is a stable fixed point of this phase (see Eqs. (f2"4|)-([26|)). 
For low temperatures (T < T c ) the system is in the ferromagnetic phase with m x ^ and, hence, 
i^l. Thus, at the critical point (T = T c ) the fixed point x c = 1 loses its stability 

u=i= i- (27) 

Here we demand that coefficients a k should be taken such that Eq.(p^) has x = 1 solution. 
For antiferromagnetics the critical temperature can be found from Eq. fl24|) and the following 
condition 

\ x=x = -i. (28) 

From definitions (|I3|)-(|T5|) for stable fixed points labeling rows (columns) of the matrixes A, A' 
and M as +|, — ^ from left to right (from up to down) we have 

M = 

{ x^Z P k = a k+1 C k p _ 2 x^-^ x^E P k = a k C k p _ 2 x^-^ 

A=(\_ ±1 and A'=( l ^ V (29) 

\ x 2 x 7 ^ / V —x 2 x 7 ^ / 



Now from Eqs. fll7| ) for two-spin correlation functions we can write 

g(n) = A 1 X n : (30) 
where A\ — 1 — m 2 (see Eq. |19]), A = A2/A1 and Ai, A2 are the eigenvalues of the matrix M 

Ax = E(a fc+2 + x 7 a fe+1 )C p fe _ 2 x^ 2 - fc ) 7 , (31) 

fc=0 

A 2 = x^ lP 5:(a fc -xa fc+1 )C p fc „ 2 x( p - 2 -^. (32) 

fc=0 

Let us now consider the behavior of g(n) in the critical region for the ferromagnetic models. 
If one restricts the interactions to pair interactions (J 2 + h), the Lee- Yang circle theorem p5H 
establishes that only at h = a phase transition can occur. However, with multisite interactions 
present (MSI+/1), phase transition may occur for h 7^ at T c 7^ and for h = only at 
T c = §. 

The T — h second order phase transition region ( surface for (MSI+J 2 + /i), line for (MSI+fo), 
and a point for ( J 2 + h) interactions) can be obtained from Eq. 



p-2 p-1 

(p -l)(q- 1) X>fc - afc+i)C?-2 = ]T a k+1 C k _ v (33) 

fc=0 fc=0 
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For example, the critical point for ( J 2 + h) interaction model on the Bethe lattice (p = 2) is [[13 

<1 



and for Husimi lattice (p = 3) is 



(J2 



{J* 



im. 

2 g-2 



1 In 25+1 



4 2g-3 

Let us now calculate the correlation length function £. Substituting x c = 1 into Eqs. (|30|)-(|32|) 
and taking into account Eq. (|33|) for correlation function at the critical region g c (n) we obtain 

g c (n) = [(p-l)(q-l)]- n . (34) 

From the definition of the dimensionality presented in Sect. 1 follows that for n 3> 1 



d 



n 



In n 



ln[(p-l)(g-l)] 



Hence, g c {n) can be written as 

g c (n) ~ 

Thus, near the critical region the correlation function has the form 



n 



where the correlation length £ is 



SW = n exp I — 



In 

(p-l)(q-l)X 



(35) 
(36) 

(37) 



(3J 



The obtained formulas (p6|)-(p8[) are in good agreement with the general behavior of correlation 
functions described in Ref. [13]. From Eqs. (pT[) , ([32]) and (|38|) follows that near the critical 
region the correlation length £ has singularity in the form 



-i 



where 



and t 



£ ~ £(Jp, J 2, h)\ t=0 t 



— j i , h , Efe=o(4+2 + 4+1 - (P - - !)K - a fc+l))Cp-2 



(39) 



J2k=l( a k+2 + ak+i)Cp_ 2 



T-Tc 



■§}CLk\t=o- Thus, the correlation length £ of ferromagnetic models on the 
Bethe-like lattices increases as the critical region is approached according to £ ~ t~ u with 
critical exponent v — 1. 

Further, let us calculate the bulk susceptibility per lattice site x — m our case, 

because we have the exact expression for correlation functions, it is easier to calculate \ using 
the fluctuation relation [13 

x = (iW) -1 (40) 
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where the summation goes over all spins, is the correlation functions of spins on sites i,j 
and Ni is the total number of sites on a lattice. Let us now consider the Eq. ( [TOD in detail 



(41) 



Taking into account the translational invariance and equivalence of all sites of Bethe-like lattice 
(see Sect. 2) we have 

(Si) = (Sj), 



N, 



Substituting this relation into Eq. (f4T|) we obtain 



X 



7^ Jim E 4 (<5bS,> - (So) 2 ) = /? lim E C^(n), 



where C n is the number of sites on the nth-shell and I is the number of shells on a lattice. From 
Eqs. (0) and fl5D| ) one can easily find 



X = P(l- m\) lim 

t— >oo 



/3(1 — m 2 ) lim 



l + g(p-l)AX)[(p-l)(9-l)A] 



n=l 



1 + A(p-1) 



l-(p-l)(g-l)A 



?(p- 1)A 



[(P-1)(?-1)A] 



l-(p-l)(g-l)A 



It is now obvious that for (p— l)(g— 1)A > 1 the susceptibility \ diverges and for (p— l)(q— 1)A < 
1 have the form 

/3(l-m 2 )(l + A(p-l)) 
X l-(p_l)( g _l)A ' 1 1 

From Eqs. ( p9|) and (P2|) for the general ferromagnetic MSI spin-^ model in an external magnetic 
field one can easily establish the relation between the susceptibility \ and the correlation length 
£ in the critical region x ~ £• 



4 Concluding Remarks 

In this paper we have applied the general Bethe-like lattice approximation to studying the 
correlation functions of the most general multisite interaction (MSI) spin-S 1 Ising models in 
an external magnetic field. Deriving the recurrence relations and using dynamical systems 
theory the thermodynamic properties of the system are studied by the investigation of the 
behavior of the fixed points of the corresponding map. Exact formulas for the correlation 
functions g(n) (spin-S" Ising models), and for the correlation length £, the bulk susceptibility 
per lattice site \ (spin-- Ising models) have been derived for arbitrary h and T. A singular 
behavior of the correlation length £ near the critical region with critical exponent v — 1 has 
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been obtained exactly. We also showed that near the critical region the magnetic susceptibility 
X is proportional to the correlation length £. 

From Eq. fl34|) one can see that in the critical region the correlation functions tend to zero 
as n — > oo and one may try to interpret this as absence of correlation on the Bethe-like lattice 
models. But this is not correct, because as we can see from Eqs. fl35|) and fl36|) the fact that at 
T = T c g(n) — > is a consequence of infinite dimensionality of the Bethe-like lattices embedded 
in the Euclidean space. As we can see from Eq. (^) correlation length in fact has singularity at 
T = T c . Hence, the second order phase transitions of models defined on the Bethe-like lattices 
have the same "attributes" as the models on regular lattices. As showed by Gujrati |TI| Bethe 
and Bethe-like lattice approximations give qualitatively correct results even when conventional 
mean-field theories fail. Furthermore, due to various options regarding the connection of the 
basic building blocks, using Bethe-like lattices one can get better approximations with respect 
to mean- field and usual Bethe lattice calculations, especially in the presence of frustrations || . 

In conclusion we want to note that the method presented here is quite general and may be 
used for other models also, for example the Potts model, the multilayer Ising model etc. 
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